CONTINUED FRACTIONS AND THE SECOND KEPLER LAW 



OLEG KARPENKOV 

Abstract. In this paper we introduce a link between geometry of ordinary continued 
fractions and trajectories of points that moves according to the second Kepler law. We 
expand geometric interpretation of ordinary continued fractions to the case of continued 
fractions with arbitrary elements. 
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Introduction 

In classical geometry of numbers the elements of an ordinary continued fraction for a 
real number a > 1 are obtained from a sail (i.e. a broken line bounding the convex hull 
of all points with integer coefficients in certain cone). In present paper we find broken 
lines generalizing sails to the case of continued fractions with arbitrary elements. This in 
its turn leads to the definition of "infinitesimal" continued fractions, whose sails would be 
differentiable curves. Such "infinitesimal" continued fractions are defined by two density 
functions: areal and angular densities. 

The areal density function has a remarkable physical meaning. Consider an observer 
at the origin and let the body move along the curve with the velocity inverse to the areal 
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density function. Then the body moves according to the second Kepler law with respect 
to the observer, i.e. with constant sector area velocity. 

This paper is organized as follows. In the first section we study the classical case of 
ordinary continued fractions. In Section 2 we expand the notion of the sail to the case 
of continued fraction with arbitrary elements. Further we show how to write continued 
fractions starting with broken lines. We generalize the proposed construction of sails to 
the case of curves, give an analog of continued fractions, and show several examples in 
Section 3. 

For a nice reference to general theory of continued fractions we suggest the book [7]. 
Several works are devoted to geometry of continued fractions (e.g. [3], [9]) and to their 
generalizations to multidimensional case ([2], [8], [4], etc). Notice that the case of broken 
lines with integer edges discussed in [5] is a particular subcase of geometric definitions 
introduced in Section 2. 

The author is grateful to V. I. Arnold for constant attention to this work, and to 
Technische Universitat Graz for hospitality and excellent working conditions. 

I. Geometry of ordinary continued fractions 

In this section we briefly introduce geometric aspects of ordinary continued fractions. 
We start with general notions of continued fractions. For arbitrary sequence of real 
numbers (a , ai, . . .) the continued fraction is an expression 
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denoted [a , a±, a 2 , . . .]. In case of a finite sequence we get some real number (or some- 
times oo). In case of infinite continued fraction the expression means a limit of a sequence 
([ao, . . . ,a n }) while n tends to infinity. Notice that such limit does not exist for all se- 
quences. A continued fraction is called ordinary if ao is integer and the rest elements 
are positive integers. A finite continued fraction is odd (even) if it contains odd (even) 
elements. 

Proposition 1.1. A rational number has a unique odd and a unique even continued 
fractions. 

An irrational number has a unique infinite continued fraction. □ 

We continue with several definitions of integer geometry. A point is said to be integer 
if all its coefficients are integers. The integer length of a segment AB is the number 
of integer points inside the segment plus one, denote it by U(AB). The integer sine of 
the angle ABC is the index of the integer sublattice generated by the integer vectors of 
the segments BA and BC in the whole lattice, we denote in by lsin(ABC). For more 
information on lattice (in particularly integer) trigonometry we refer to [5] and to [6]. 

Let C be a cone with vertex at the origin. Take the convex hull of all integer points 
except the origin in C. The boundary of the described convex hull is a broken line together 
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a = UiAoAt) = 1; 

ai = lsin/A AiA 2 = 2; 

a 2 = U^Az) = 2. 

7/5 = [1,2,2]. 



Figure 1. The broken line A A 1 A 2 is the sail for C 7 / 5 . 

with one or two rays in case if there are some integer points on the edges of the cone. The 
broken line in the boundary of the convex hull is called the sail of C . (In some literature 
the sail is the whole boundary of the convex hull, for us it is more convenient to exclude 
the rays from the definition of a sail.) 

Consider a positive real number a. Denote by C a the cone with vertex at the origin 
and edges {(t,0)\t > 0} and {(t,at)\t > 0}. The sail for C a is a finite broken line if a is 
rational and one-side infinite broken line if a is irrational. 

So let the sail be a broken line A Q Ai . . . A n (A AiA 2 . . .). Denote 

a2fc-i = W(A k A k+1 ), 
a 2k = lsin(A fc _iA fe A fc+ i) 

for all admissible k. The sequence (ao, a±, . . . , a 2n ) (or (a , a±, a 2 , ■ ■ .)) is called the lattice 
length-sine sequence for the cone C a (or LLS-sequence for short). 

The connection of geometric and analytic properties of LLS-sequence is introduced by 
the following theorem. 

Theorem 1.2. Let a > 1 and (a , ai, . . . , a 2n ) (or (a , a-y, a 2 , . . .)) be the LLS-sequence 
for C a . Then 

a = [ao, ai, . . . , a 2n ] (or respectively a = [ao, a±, ...]). 

Proof. This theorem is a reformulation of a Proposition 1.7. a from [5] for finite continued 
fractions and Theorem 2.7.a from [6] for the infinite continued fractions. Finite case is 
also a particular case of Corollary 2.6. So we skip the proof here. □ 

2. Continued fractions with arbitrary coefficients 

In this section we generalize geometry of ordinary continued fractions to the case of 
continued fractions with arbitrary elements. We show a relation between odd or infinite 
continued fractions and broken lines in the plane having a selected point (say, the origin). 
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PQ = ct2kOA k ; 
A k A k+1 = a 2k+ iA k Q. 



Figure 2. Construction of A k+1 . 

We conclude this section with a few words about conditions for a broken line to be closed 
in terms of elements of the corresponding continued fraction. 

Further we use the following notation. For a couple of vectors v and w denote by |i> x w\ 
the oriented volume of the parallelogram spanned by the vectors v and w. 

2.1. Construction of broken lines from the elements of continued fractions. In 

this subsection we give a natural geometric interpretation of an odd or infinite continued 
fraction with arbitrary elements. It would be a broken line defined by the positions of 
the first vertex and the selected point O, direction of the first edge, and the continued 
fraction. 

So consider a continued fraction [a , . . . , a 2n }. We are also given by the vertex A , 
selected point O, and the direction v of the first edge. We construct all the rest vertices 
A k inductively in k. 

Base of induction. For the second vertex we take 

A 1 = A + Xv, 
where A is defined from the equation \OA x OAi\ = a±. 

Step of induction. Suppose now we have the points A , . . . , A k , for k > 1, Let us get 
A k+1 . Consider a point 

P = A k + ^±^A k _ lAk . 

a 2k-l 

In other words P is a point in the line A k _iA k such that the area OA k P equals 1. Let 

Q = P + a 2k OA k . 

Finally the point A k+ i is defined as follows (see on Figure 2) 

A k+ i = A k + a 2k +iA k Q. 

Let us now explain a geometric meaning of the elements of continued fractions in terms 
of characteristics of the corresponding broken line. 

Proposition 2.1. The following holds 

a 2k+ i = \OA k x OA k+1 \, k = 0,...,n, 
\A k A k _i x 

a 2k = , fc = l,...,n. 

a2fc-l02A:+l 
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Proof. We prove this statement by induction in k. 

Base of induction. From the definition of the point A 1 we get 

\OA x OAi\ = ai. 

Step of induction. Let the statement be true for k—1, we prove it for k. 
First, we verify the formula for a 2k +i' 

\OA k x OA k+1 \ = a 2k +i\OA k x OQ\ = a 2k+ i\OA k x OP\ = -^±i|OA fe _i x OA k \ = a 2k+1 . 



Q-2k-l 



The last equality holds by induction. 

Second, for a 2k we have 

L4fcAfc_i x A k A k+1 \ = \A k A k ^ x A k Q\ 

a 2k-l a 2k+l a 2k-l 

\OA k _i x OA k \ = a 2k . 



\PA k x PQ\ = a 2k \OA k x OP\ 



a 2k-l 

The step of induction is completed. □ 

Example 2.2. Let us construct a broken line having the first vector Aq = (1,0), the 
direction v = (1, 0), and the continued fraction [a, b, c]. Then we have 

A 1 = (l,a). 

Further we find the corresponding points P and Q: 

P = (l,l + a), Q = (l + b, 1 + a + ab). 

Finally we get 

A 2 = (1 + be, a + c + a&c). 

2.2. Inverse problem. Now suppose we have a point O and a broken line A . . . A n such 
that for any k the points O, A k , and are not in a line. Let us extend the definition 
of the LLS-sequence for this data. 

We use equalities of Proposition 2.1 to define the elements: 

a 2k+1 = \OA k x OA k+1 \, k = 0,...,n; 

„ _ \A k A k _i x A k A k+ i\ , 

°2fc- a 2fc _ia 2fc+ i ' fc-l,...,n. 

We call the sequence (a , . . . , a 2n ) the LLS-sequence of the broken line with respect to the 
point O, and [ao, • • • , a 2n ] — the corresponding continued fraction. 

Proposition 2.3. Let A . . . A n and B . . . B n be two broken lines with LLS-sequences 
(a , . . . , a 2n ) and (b , . . . , b 2n ) respectively. Suppose the first broken line is taken to the 
second by some operator in SX(2,R) with determinant equals A. Then we have: 

a 2k = Xb 2k , k = 1, ... ,n 
fl2fc+i = jb 2k+1 , k = 0, . . . , n 

Proof. The volume of any parallelogram is multiplied by A, then the statement follows 
directly from formulas of Proposition 2.1. □ 
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2.3. On geometric meaning of corresponding continued fractions. For this sub- 
section we fix the point O to be at the origin. 

Consider a continued fraction [ao, Oi, . . . , a k ] as a rational function in variables do, . . . , a&. 
This rational function is a ration of two polynomials with non-negative integer coefficients, 
denote them by P k and Q k . Actually the polynomials P k and Qk are uniquely defined by 
the condition 

Pk\ a 0i °1) • • • j a k) r n 

— — r — [ao, ai, . . . , a k \ 

LJ k {a , a>i, ■ ■ ■ , OfcJ 

and the condition that the coefficients of both polynomials are non-negative integer coef- 
ficients. 

Remark 2.4. Notice that the last condition is equivalent to the condition that the poly- 
nomial -Pfc(a , ai, . . . , a k ) contains a monomial a a>i • . . . • a k with unit coefficient. 

Theorem 2.5. Let A . . . A n be a broken line such that A = (1, 0), and A\ = (1, ao) is 

collinear to the vector (0, 1). Suppose its LLS-sequence is (ao, a±, . . . , a,2 n )- Then 

A n = (<?2n+l( a 0? 0,1, ... , a 2n ), -P2n+l( a 5 a l; • • • 5 a 2n)) ■ 

Proof. We prove this statement by induction in n. 

Base of induction. If the broken line is a segment A Ai with LLS-sequence (ao) then 
A 1 = (l,a ). 

Step of induction. Suppose the statement holds for all broken lines with k vertices, let 
us prove it for an arbitrary broken line with k + 1 vertex. 

Consider a broken line Aq. . .A^ with LLS-sequence (ao, . . . , a,2k)- Let us apply a linear 
transformation with unit determinant taking A 1 to (1,0) and the line A 2 Ai to the line 
x — 1. This transformation is uniquely defined by all these conditions, it is 

aoai + 1 — ai 
-a 1 

Denote the resulting broken line by B P>i . . . B k . By Proposition 2.3 all the elements of 
the LLS-sequence for B Bi . . . B k are the same. By the assumption of induction we have 

B k = (Q2k-l(^2, • • • ? 0-2k)i P2k-l{0'2, ■ ■ ■ , 0,2k))- 

Denote the coordinates of B k by q and p respectively. Then we have 

A k = r _1 (-Bfe) = (p + aiq,a p + (a ai + l)q). 

The polynomials satisfy 

QqP + (oqQi + l)g _ a 1 _ P2k+i(ao,ai, ■ ■ ■ ,a>2k) 

P + a x q Q i + f Q2k+i(a ,a l , . . . ,a 2k ) 

Notice that the polynomial aop + (aoOi + l)q has a unit coefficient in the monomial 
aoai • . . . • a,2k coming from aoaig. Therefore (see Remark 2.4), a^p+ (a^ai + l)g coincides 
with P2fc+i(oo, ai, . . . , a 2 fc) and (p + a±q) coincides with Q2k+i(o , a±, . . . , a 2k )- So we are 
done with the step of induction. This concludes the proof of the theorem. □ 
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In particular we get the following corollary. In the classical case it forms the basis of 
geometry of ordinary continued fractions. 

Corollary 2.6. Let A ...A n be a broken line such that A = (1,0), and A = (l,ao). 
Suppose that the corresponding continued fraction is a = [a , ai, . . . , 02 n ] and A n = (x, y) . 
Then 

y 

— = a. 

x 

(If the corresponding continued fraction has an infinite value, then x/y = 0.) □ 

This corollary implies the following statement. 

Corollary 2.7. Let A ...A n and B . . . B m be two broken lines with B = A , such 
that the vector of the first edges either have the same direction if ao/^o > or opposite 
otherwise. Suppose the corresponding continued fractions coincide: 

[do, • • • , a2n] = [bo, • • • , &2m]- 

Then the points A n , B m , and the origin O are in a line. 

Proof. Consider an SL(2, M)-operator taking Aq to (1,0) and A\ to the line x — 1. By 
Proposition 2.3 the continued fractions for both broken lines are not changed. Hence by 
Corollary 2.5 the points A n , B m , and the origin are in a line. □ 

Remark 2.8. (On closed curves.) How to find that a certain continued fraction defines a 
closed curve? From Theorem 2.5 we see that a broken line defined by an LLS-sequence 
(ao, ai, . . . , a2n) with A = (1, 0) and A Ai being collinear to (0, 1) is closed if and only if 

Q2n+i(ao,a u . . . ,a 2n ) = 1 and P 2n+ i(a , a 1: . . . , a 2n ) = 0. 

So, these two polynomial conditions on the elements of the LLS-sequence are necessary 
and sufficient conditions for the broken line to be closed. 
Notice that the condition P n+ i = can be rewritten in the following nice form 

[do, ai, . . . , a 2n ] = 0. 

The condition P n+ i = was introduced in [5] for certain broken-lines with integer vertices. 

Example 2.9. Let us study an example of broken lines consisting of three edges. These 
curves are defined by continued fractions of type [00,01,02,03,04]. Then the conditions 
for a broken line to form a triangle are as follows: 

J 0001020304 + a aia 2 + a aiO4 + 000304 + 020304 + 00 + 02 + 04 = 
1 01020304 + Oi02 + Oi0 4 + 0304 + 1 = 1' 

(See on Figure 3.) 

There is one problem which is interesting in the frames of this section. Suppose we 
have a broken line and two distinct points 0\ and 2 . Then we have two LLS-sequences 
for the same curve with respect to 0\ and 2 . Study the conditions on the initial data 
(i.e., LLS-sequences, positions of the first points of the broken lines, and direction of the 
first vector) that define congruent broken lines. 
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A 3 A 1 -2 




[2, -1, 3, -2, 1] = [1, -2, 2, -1/2, -4] = ^ 

Figure 3. Examples of broken lines and their continued fractions. 



3. DlFFERENTIABLE CURVES 

Now let us study what happens if we consider a curve as a broken line with infmitesi- 
mally small segments. It turns out that the LLS-sequence "splits" to a couple of functions 
which we call areal and angular densities. We introduce the necessary notions and discuss 
basic properties of these functions. In particular we show that the areal density is inverse 
to a velocity of a point defined by the second Kepler law. 

In this section we suppose that the curves has a natural (unit length) parametrization. 

3.1. Definition of areal and angular densities. Consider a curve 7 of class C 2 with 
an arc-length parameter f. Let us define the areal and the angular elements at a point 
similar to the discrete case. 

Definition 3.1. The areal density and the angular density at t are respectively 

|0 7 (f)x0 7 (f + e)| 



A(t) 



lim 



and 



B(t) 



lim 

e- 



= |07(f) x 7(f) I 



| 7 (f) 7 (f-£) X 7 (f) 7 (f + £)| 



-0 e|Oy(t - e) x 07(f) 1 1 7 (f) x 7 (f + e) | ' 
Let us give geometric interpretations for the functions A and B. We start with A. 

Proposition 3.2. (Relation with the second Kepler law.) Suppose that a body 
moves by a trajectory of a curve 7 with velocity 1/ A. Then the sector area velocity of a 
body is constant and equals 1. 



Proof. The proof follows directly from the definition. 



□ 



Instead of giving a geometrical interpretation of B, we prove the following formula for 
A 2 B. For a given curve 7 denote by n(t) the signed curvature at point f. 

Proposition 3.3. Consider a point 7(f) of a curve 7. Let the vectors 07(f) and 7(f) be 
non-collinear. Then the following holds. 

A 2 (t)B(t) = K (t). 
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Proof. We have the following 

|07(g x Q 7 (t + e) I \ 2 | 7 (t) 7 (t - e) x 7 (t) 7 (t + e)| \ 

e ' ; e\0 1 (t-e)x0 1 (t)\\0 1 (t)x 0^ + 8)1 J 

-lim l7^7(t-g) x 7 (t) 7 (t + £ )| 
e^o e 3 

Notice that 

\>y(t)>y(t-e)\=e + o(e), 
| 7 (t)7(* + e)l =e + o(e), 
sin( 7 (t - e) 7 (t)7(* + £ )) = + 

Therefore, for the volume of the corresponding parallelogram we get 

| 7 (t)7(* - e) x t(*)t(* + e)l = + 
Hence, A 2 (t)B(i) = «(*). □ 

Now we prove the theorem on finite reconstruction of a curve (i.e., in some small 
neighborhood) knowing the areal density and a starting point. This is analogous to the 
algorithm that finds a broken line by the elements of the corresponding continued fraction 
described in Subsection 2.1. The significant difference to the discrete case is that we do 
not need to know the angular distribution function. 

Theorem 3.4. Suppose we are given by the points O and 7 (to) an d the areal density 
A(t ). 

- If \A(t )\ > |O 7 (t )|; then there is no finite curve with the given data. 

- If |O 7 (t )| > |^4(to) I > 0, then the curve is uniquely defined in some neighborhood 
of the point j(to). 

Remark 3.5. Notice that A 2 (t)B(t) defines the oriented curvature. Therefore, if one knows 
the functions A and B then the curve is uniquely reconstructed until the time t where 
the vectors O 7 (t ) and 7 (to) are collinear, or in other words where |A(i )| = 0. 

Proof. Consider a system of polar coordinates (r, </?) with the origin at the point O. To 
get the curve we should solve the system of differential equations: 

f r 2 ip = A 

\ r 2 + r 2 ip 2 = 1 ' 

This system is equivalent to the union of the following two systems: 

By the main theorem of theory of ordinary differential equations (see for instance in [1]) 
this system has a finite solution if \r\ > \A\ > 0. This concludes the proof. □ 



A 2 (t)B(t) = lim 
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Let us say a few words about density functions and their broken line approximations. 
Let 7(t) be a curve with arclength parameter t G [0,T] and densities A(t) and B(t). 
For an integer n consider a broken line j n = Ao, n • • • A n ^ n such that Ai tTl = 7(-T). Let 
the corresponding LLS-sequence be (a , n , . . . , a-2n,n)- Denote by A n and B n the following 
functions 

A n (t) = a 2 \nt/T\+i,n, and B n {t) = a 2 \nt/T\,n- 

Theorem 3.6. Let 7 be in C 2 . Then the sequences of functions (A n ) and (B n ) pointwise 
converge to the functions A and B respectively. 

Proof. This follows directly from the definition of density functions and Proposition 2.1. 

□ 

It is interesting to investigate the inverse problem, it is still open now, we formulate it 
in the last subsection. 

3.2. Example of curves and their continued fractions. In this subsection we calcu- 
late the areal and angular densities for straight lines, ellipses, and logarithmic spirals. 

Example 3.7. Lines. Let us study the case of lines. Without lose of generality we 
consider the point O to be at the origin and take the line x — a. Then the corresponding 
densities are 

A(t) = a and B(t) = 0. 

Example 3.8. Ellipses and their centers. Consider an ellipse ^ + |j- = 1 with 
a > b > 0. Let O be at the symmetry center of the ellipse i.e. at the origin. Then the 
areal and angular densities are as follows 

Ait) = == ah === and B(t) = 



Va 2 sin 2 1 + b 2 cos 2 t abV a 2 sin 2 1 + b 2 cos 2 1 

Notice that here we get the constant function for the ratio: 

A (t) 0,2 



B(t) 



= a 2 b 2 . 



Example 3.9. Ellipses and their foci. As in the previous example we consider an 

2 2 

ellipse ^2 + 52 = 1 with a > b > 0. Let now O be at one of the foci, for instance at 
(— \/a 2 — b 2 , 0). Then the densities are as follows 

. . . ab + b\J a 2 — b 2 cos t , „ . . a 

A(t) = and B(t) = ^ 2 . 

V a 2 sin 2 1 + b 2 cos 2 vt bV a 2 sin 2 1 + b 2 cos 2 t(a + cos t\J a 2 - b 2 ) 

Remark on the Kepler planetary motion. If we put the Sun at the chosen focus and 
consider a planet whose orbit is the ellipse, then according to three Kepler laws the planet 
will move with velocity \/A{t) at any t. Here the constant A is defined from the third 
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Kepler law: the square of the orbital period of a planet is directly proportional to the cube 
of the semi-major axis of its orbit, or in other words 

e 

a 3 a\ ' 

where T is the period for our orbit, and T e and a e are respectively the period and the 
semi-major axis for the Earth. Denote by L the length of the ellipse, i.e., 

tt/2 

b 2 ' 

— ) cos 2 t dt 
a 2 




Since T — |A| J* \1/A(t)\dt, we get 

A = : 



, 3 

71 / a \ 3 



Jo |lM(t)|^t V«e 

We skip a description for parabolas and hyperbolas, they are similar to the case of 
ellipses. 

Example 3.10. Logarithmic spirals. Consider a logarithmic spiral 

{(ae bt cost, ae bt sin t)\t G R}. 
Then the densities for this spiral are as follows 



A(t) = . and B{t) 



^/WTl a 3 
It is interesting to notice that for the spirals we have 

A\t)B{t) 1 



62 + 1' 

i.e., the products are constant functions. 

Notice that if A 2 B is a constant function, then the curvature is constant, and hence we 
get circles. What do we have if AB (or A) is constant? 

3.3. Open problems. We conclude this section with two open problems concerning the 
density functions. We start with a question on convergency that is in some sense the 
inverse problem to Theorem 3.6. 

Problem 1. What properties should have the LLS-sequences of broken lines if their 
sequence converges to certain curve. 

The second problem comes from Remark 2.8 on closed broken lines. 

Problem 2. What are the conditions on the functions A{t) and B(t) for the resulting 
curve 7 to be closed? 
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